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We study the flux-driven superconductor-metal transition in ultrasmall cylinders observed exper- 
imentally by Liu ei.a/.pjj. Where T c — > 0, there is a quantum critical point, and a large fluctuation 
conductivity is observed in the proximate metallic phase over a wide range of T and flux . However, 
we find that the predicted (Gaussian) fluctuation conductivity in the neighborhood of the quantum 
critical point is 4 orders of magnitude smaller than observed experimentally. We argue that the 
breakdown of Anderson's theorem at any non-integer flux leads to a broad fluctuation region re- 
flecting the existence of "rare regions" with local superconducting order. We calculate the leading 
order correction to the conductivity within a simple model of statistically induced Josephson coupled 
local ordered regions that rationalizes the existing data. 



The physics of the zero temperature "superconductor- 
insulator" transition in low dimensional systems is im- 
portant in its own right, and as a paradigm for quantum 
phase transitions in fermionic systems 2] . It is therefore 
embarrassing that it remains an unsolved and contraver- 
tial problem. Most studies have focussed on cases in 
which the transition is either driven by disorder or a com- 
bination of disorder and magnetic field. However, these 
important studies often reveal the existence of anomalous 
metallic (dissipative) phases where superconductivity is 
destroyed, rather than the expected insulating phases. 
(Some representative reviews of this subject are con- 
tained in Refs. HSU Ho 

Moreover, it is presently unclear whether the range of 
parameters (e.g. temperature, T, and magnetic field) 
over which quantum critical fluctuations dominate the 
behavior of such systems is large or small - this later 
uncertainty likely has more general implications for the 
issue of whether or not quantum critical fluctuations un- 
derly a host of "non-Fermi liquid" properties seen in 
highly correlated materials from heavy fermion metals 
to the high temperature superconductors. (See, for ex- 
ample, Refs. HQ, 00.) 

We study here a concrete and particularly simple ver- 
sion of this problem that was recently studied experi- 
mentally by Liu et. al. lj : As shown in the inset of Fig. 1 , 
the system consists of a thin insulating cylinder coated 
with a simple BCS metal. Magnetic flux of magnitude 4> 
is threaded along the central axis. The quantum phase 
transition is driven by the balance between the kinetic 
energy stored in the flux induced diamagnetic currents 
and the condensation energy, i.e. this represents the ex- 
treme limit of the Little-Parks experiment. Experimen- 
tally, it was observed that at a critical value of <f> = 4> c 
there is a continuous quantum phase transition from a su- 
perconductor to a metallic phase, where 4> c < &o/2 and 
$o = hc/2e is the superconducting flux quantum. How- 
ever, even well beyond the critical flux, where <j> ~ $0/2, 
the T = resistance is ~ 0.3 of its normal state value, 
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FIG. 1: Schematic phase diagram for the flux induced quan- 
tum phase transition in a thin walled superconducting cylin- 
der. The superconductivity survives inside each dome, while 
at temperatures above the mean-field T c , the system is a 
Fermi liquid. The "Lifshitz phase" at intermediate <f> is dom- 
inated by superconducting fluctuations in rare regions with a 
larger than average local T c . The dotted line denotes half the 
flux quantum $o/2 = hc/4e. Inset - sketch of the experimen- 
tal geometry. 



indicating the persistence of substantial superconducting 
fluctuations. 

We have found that the value of (j> c can be obtained 
with surprising accuracy from simple mean-field theory. 
However, standard^] small amplitude (Gaussian) fluc- 
tuations about the mean-field solution are at least four 
orders of magnitude too small to account for the observed 
fluctuation conductivity. Rather, we argue that a break- 
down of Anderson's theorem (due to the breaking of time- 
reversal symmetry by the applied flux) combined with 
the physics of rare regions characteristic of disordered 
systems (akin to the physics of Lifshitz tails in semicon- 
ductors) provide a natural explanation for at least the 
bulk of the observed fluctuation effects. (The nature of 
the non-Fermi liquid groundstate at T = remains an 
intriguing, but still unsolved problem.) 



For </> = 0, the system is time reversal invariant, hence 
Anderson's theorem implies that T c (£,(f> = 0) is inde- 
pendent of the mean-free path, £, over a broad range of 
I > kp , where Hf is the Fermi wave number. Hence, 
although there are always mesoscopic variations in the 
local values of £, T c remains uniform, and indeed the ob- 
served resistive transition is very sharp. 

On the other hand, at non-zero <fi, time reversal is no 
longer a symmetry, and hence T c (£, cf) varies with £; cor- 
respondingly, at T = 0, the critical value of the flux is (as 
shown in Eq. a strong function of £ (For related issues 
see Refs. 0,EIE3)- As a consequence, even when 4> is 
greater than the average value of (j> c , local statistical fluc- 
tuations of £ lead to rare regions with a non-zero local T c . 
At not too low temperatures, the effect of these regions 
on the metallic conductivity can be calculated in a per- 
turbative expansion in the Josephson coupling between 
them, Ej (See Eq. ^3) The phase dynamics are dom- 
inated by their coupling to the dissipative environment, 
and as a result they can make a substantial contribution 
to the conductivity even when T > Ej. The result is a 
broad region of the <p — T plane, which we have called the 
"Lifshitz-phase," in which the resistance is reduced from 
its value in the Fermi liquid phase. 

The Hamiltonian for a single electron confined to the 
cylinder or radius Rq is 
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where we work in a gauge in which A points in the az- 
imuthal direction with Ag = <f>/(2-Kr) w 4>/(2irR ). Since 
most of the flux threads through the insulating midsec- 
tion, the magnetic field is small despite the large value 
of the flux, so the Zeeman coupling is negligible. As long 
as the magnetic length is longer that the thickness of the 
thin film coating, any additional orbital effects associated 
with the curvature of the electronic paths can be ignored 
as well. 

The disorder averaged pairing susceptibility of a de- 
generate assembly of such electrons is given by 

Xp (k,ztt) =N(0)k B Tj2 I • fwilP 
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where q = ( 1 - k, ( = (^ /<j>)Ro, and g(q, A) w 
— % -i D ni 4 for diffusive electrons with the diffusion con- 
stant D = Vpl/3. The above integrals can be evaluated, 
yielding 
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As usual, in the weak coupling limit the condition 1 = 
5Xp(0,0) determines the critical temperature. Near the 
quantum critical point, Xp(k, iO) can be expanded in 



powers of T 2 and we get 
Xp (k,m) = ^N(0) (in 
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which means that near the quantum critical point 
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gives the width of the superconducting domes in Fig. 
1 centered around each integer flux <f> — n$ j so long 
as <j) c < $o/2; otherwise, T c never vanishes, but rather 
reaches a minimum, non-zero value for <fi = §q/2. 

For Rq — 60nm, the superconducting coherence length 
in clean Al £o = 1600nm, and the mean-free path, derived 
from the normal state conductivity, £ — 12nm, Eq. (@J 
yeilds a theoretical estimate, <fi c = 0.42$o, which agrees 
with experiment |H to better than 5% accuracy. The 
agreement with the experiment, in addition to the nar- 
rowness of the transition regime in the absence of flux 
, demonstrate that the superconductivity, as assumed, 
has a bulk character and the finite size of the cylinder 
radius only serves to determine <p c . 

Inspired by the success of the mean field theory for 
4> c , we now turn to the question of the character of the 
disordered state when <&o/2 > 4> > 4> c - In particular, we 
calculate the extra conductivity on the disordered side. 
The action governing the dynamics of the fluctuating or- 
der parameter is 

S = J [A- r _ r ,(r-r')At(r,r)A(r',r') + foo|A(T,r)| 4 ] 

(5) 

where A is the charge 2e pairing field, and the integral 
ranges over space-time; K(uj,q) = g^ 1 — x p (w,q) and 
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where ip( n ~ lS> is the n— th logarithmic derivative of the 
Gamma function. 

Within the Aslamazov-Larkin theory fluctuation con- 
ductivity in d-dimensions is 
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where the mass m is given in Eq. @. This contribution 
vanishes at T — except possibly at the quantum critical 



point. Additional effects coming from Maki-Thompson 
and density of states diagrams were computed by Lopatin 
et.al. 15]. In addition, Ref.flij included frequency de- 
pendence of the boson vertex which gives a contribution 
at T=0, but contrary to the experiment Q leads to a 
weak negative magnetoresistance. Overall, the critical 
regime where the Gaussian fluctuations contribute ap- 
preciably is at least 4-orders of magnitude smaller than 
observed experimentally Q. 

In fact, in the absence of the flux, the superconducting 
transition is extremely narrow, consistent with the small 
value of the Ginzburg parameter of a low T c material, 
while near the quantum critical point the transition is 
very broad. As mentioned previously, at non-zero 0, T c 
is a strong function of both I and the cylinder radius. 
Since at the mean field quantum phase transition, the 
Harris criterion is not satisfied, i.e. vd < 2, one expects 
a disorder driven instability of the clean critical point. 



Thus, the ensemble fluctuations of £ and Rq in the system 
lead to randomness in the mass of the effective quantum 
Ginzburg Landau theory. As a consequence, rare regions 
in space where 4> < 4> c develop local superconducting 
order and can contribute anomalously to the "normal" 
state conductivity. 

To be explicit, we model the disorder of the critical 
coupling constant by a set of Ni mp delta function "impu- 
rities" of strength — \u\ positioned at random throughout 
the sample. This model, while highly simplified, allows 
analytic treatment in the limit of small impurity con- 
centration, i.e. the regime marked by the onset of the 
anomalous contribution to the conductivity. 

Since w, the thickness of the Al thin film coating, as 
well as the cylinder radius, Rq 7 are much less than the 
coherence length, the order parameter can be taken to 
be a function of the position along the cylinder, z, only. 
The resulting Ginzburg-Landau functional is then 
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where £ = ($o/4>)Rq, the Fourier transform of the inverse 
propagator G _1 (w) = <^ 2 \oj\/D, the contact repulsion b = 
bo/N(0), and the average "mass" 
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The saddle point of the static action with a single "im- 
purity" at the origin is 
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where A = tanh -1 [2m£/|u|]. The order parameter takes 
its maximum value, Ao = \J (u/2£) 2 — m 2 / s/b, at the 
origin. Note, that this solution minimizes the action only 
for 2m( < \u\, where A is real. 

If the concentration of such impurities, Ni rnp /L <C 
m/C, (L is the length of the cylinder), is small we can 
adopt as an ansatz for the saddle-point solution 

A(t,z) = J2 e ie ^A lmp (z- Zj ). (11) 
i=i 

Substituting the above equation into a Bogoliubov- 
de Gennes equation and integrating out the fermions 
to leading order in the Josephson coupling strength, 
Ej, gives an effective action for the phase differences 
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where a — Rq/Rn, the quantum of resistance Rq = 
h/4e 2 = 6.5k£l, So is the part of the action independent 
of the phase configuration 0, and Ej(j) is the Joseph- 
son coupling between the sites. In the dilute impurity 
limit and at finite temperature, Ej{j) is given approxi- 
mately by 



2\u\ A . 

Ejx 1 ) = — A imp (z i+1 

4 e c 



(13) 



where e" 1 = 27ri? ^CA r (0), i.e. the Josephson coupling 
falls off exponentially with m\zi+\ — Zi\jC,. 

So long as T is not too low, the correction to resis- 
tivity can be computed in powers of Ej. The most di- 
rect way is to use the Keldysh technique |l7l. with the 
appropriate Keldysh propagator matrix |l8|, G^(uj,r), 
G 12 (u,r) = G 21 (-u,r) - G R (u;,r^ G u (w,r) = 
coth [tf/2T] (Gi2(w) — G2i(w)), and the retarded prop- 
agator is G (o),r) = l/(iauj), which is regularized in 
the UV by introducing a cuttoff Tjjv ~ ^/(Abt). For a 
given disorder configuration, the resulting correction to 



the normal state resistivity is 
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where g(a) = cos[£] {^f /a T[\ - ±]r[±]/(2^a). In 
the thermodynamic limit, the arithmetic mean of E 2 (j) 
equals its (disorder) ensemble average, 

W—Y, E JW= [^1 P(l)Ej[l}, (15) 

1* imp . Jo 



where P(l) is the distribution of the junction lengths and 
Ej[l] is the Josephson coupling across a junction of length 
I. In the model of disorder that we used, P(l) is the prob- 
ability density of finding a variable with a value I among 
Nimp independent variables, each evenly distributed be- 
tween [0, L], and whose sum is constrained to be L. In 
the limit of large L, with the density rii mp = Ni mp /L 
held fixed, P(l) can be shown to be given by the prob- 
ability distribution for the occurence of the first success 
in a Poisson process, i.e. P(l) = n imp e~ nimpl . Thus, 
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So we arrive at one of the main results of our analysis (to 
+0(Ej)) 
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In a good metal, a 1 and the above expression has 
only a weak dependence on the high energy cuttoff, Tjjy 
The main point is the explicit functional dependence on 
the scale of the local ordering Ao and the density of the 
locally ordered regions expressed through the dimension- 



less parameter c 
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Note that in general both Ao and c have implicit tem- 
perature dependence. If the deviations Sm(z) from the 
mean mass m are S— function correlated gaussian random 
variables, (Sm(x)5m(y)) = ^W5(x— y), then for m ^> W, 
n, the density of Ginzburg-Landau "impurities" (see Eq. 
|SJ) is given by the cumulative density of states in the Lif- 
shitz tail below the energy m 2 . In 1 dimension this is 
given by n(m 2 ) = ±f e -ie(m 3 /3W H) Similarly, the 
strength of the impurity it, is chosen by matching the 
average negative eigenvalue in the Lifshitz tail below m 2 
with the eigenenergy of the impurity bound state. Thus 
(u/(2C)) 2 = m 2 + W/{8m() which gives A| = W/(8m£b). 

Fitting Eq. (jTo)l to the data of Ref. [l[ yields a qual- 
itatively good fit for physically reasonable parameters. 
Hence our theory provides a natural explanation for the 
broadness of the transition observed in the experiment. 
Quantitative details will require further investigation. 



We conclude with a discussion of the validity of our ap- 
proach. While the dynamics of the local ordered regions 
can almost certainly be described by the Eq. (fT2")) at high 
temperature, as the temperature is lowered, the action 
(|12fl must be modified. There are several reasons: first 
as the temperature is lowered, the locally ordered regions 
induce suppression of the single particle density of states 
at the Fermi level, thereby reducing dissipation of the col- 
lective phase variable and suppressing phase ordering. In 
addition, the simple harmonic action describing dissipa- 
tion is suspect at low temperature. Finally, the quantum 
mechanical tunneling of the order parameter amplitude 
can lead to a destruction of the phase orde ring; at low 
temperature, giving rise to a dissipative state 20]. 
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